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We analyze a one dimensional version of a model of morphogen transport, a biological process 
governing cell differentiation. The model was proposed by Hufnagel et al. to describe the 
forming of morphogen gradient in the wing imaginal disc of the fruit fly. In mathematical terms 
the model is a system of reaction-diffusion equations which consists of two parabolic PDE's and 
three ODE's. The source of ligands is modelled by a Dirac Delta. Using semigroup approach 
and L\ techniques we prove that the system is well-posed and possesses a unique steady state. 
t-h All results are proved without imposing any artificial restrictions on the range of parameters. 

AMS classification 35B40, 35092. 
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1 — 1 1 Introduction 

> 

^j- Morphogen transport (MT) is a biological process occurring in the tissues of living organisms. It 
^ is known that certain proteins (ligands) act as the morphogen - a conceptually denned substance 
t^j- which is responsible for the development and differentiation of cells. As it is proposed in the 'French 
flag model' by Wolpert [TB], morphogen molecules spread from a spatially localized source through 
O the tissue and after some time form stable concentration gradient. According to the concept of 
,_i positional signalling, receptors located on the surface of the cells, detect information about local levels 
of morphogen concentration. This information is transmitted to nucleus and cause gene activation 
which finally leads to the synthesis of suitable proteins and cell differentiation. Although the role of 
morphogen gradient in gene expression seems to be widely accepted, the exact kinetic mechanism of 
its formation is still not known (see [3] , [5] , [H] ) • 

Recently various models of MT, consisting of PDE-ODE systems, were proposed (see [TO], [I], jl]) and 
analyzed (see [6] , [7] , [15] , [13] , [H] ) . These models assume that movement of morphogen molecules 
occur by different types of diffusion or by chemotaxis in the extracellular medium. Reactions with 
receptors (reversible binding, transcytosis) and various possibilities of degradation and internalization 
(of morphogens, receptors, morphogen-receptor complexes) are also being considered. 

In [1] Hufnagel et al. proposed a model, which we denote [HKCS], of the formation of morphogen 
Wingless (Wg) in the Drosophila Melanogaster wing imaginal disc. Apart from mentioned before 
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processes, the model takes into account interaction of Wg with glypican Dally-like (Dip) - protein 
which, similarly to receptor, interacts with morphogen through association-dissociation mechanism. 
Dip molecules may also transmit Wg to each other causing the movement of morphogen particles on 
the surface of the wing disc. The interesting issue is the presence of a singular term (a Dirac Delta), 
to model the secretion of morphogen from a narrow part of the tissue which is represented in the 
model by a point source. In the present paper we analyze a one dimensional simplification of the 
[HKCS] model, which was also introduced in pi]. We intend to extend our analysis to the case of 
higher dimensional domains in the forthcoming article(s). 



1.1 The [HKCS]. ID model. 

In this section we present a one dimensional simplification of the model [HKCS] introduced in [I]. 

For L > 0, oo > T > 0, denote I L = (— L, L), dI L = {— L, L}, 1% = (0,T) x I L , (dI L ) T = 
(0, T) x dI L , I = I 1 . The following system is of our interest: 

[HKCS]. ID 

d t W = Dd 2 xx W - 7W - [kGW - k'W*} - [k R RW - k' R R*} + s6, (t, x) G i£ (la) 

d t W* = D*d 2 xx W* - ^*W* + [kGW - k'W*] - [k Rg RW* - k' Rg R* g ], (t, x) G (lb) 

d t R = -[k R RW - k' R R*] - [k Rg RW* - k' Rg R* g ] -aR + T, (t, x) G (lc) 

d t R* = [k R RW - k' R R*] - a*R*, (t, x) G l£ (Id) 

d t R g = [k Rg RW* - k' Rg R* g ] - a*R* gl (f , x) G /£ (le) 

d x W = d x W* = 0, (t, x) G {dl^ (If) 

W(0) = W , W*(0) = W*, R*(0) = R* , R* g (0) = R* g0 , R(0) = R , x G I L (lg) 

In Q W,G,R denote concentrations of free morphogens Wg, free glypicans Dip and free receptors, 
W*,R* denote concentrations of morphogen-glypican and morphogen-receptor complexes, 
R* denotes concentration of morphogen-glypican-receptor complexes. 
The model takes into account association-dissociation mechanism of 

• W and G with rates k, k' (kGW - k'W*), 

• W and R with rates k R , k' R (k R RW - k' R R*), 

• W* and R with rates k Rg , k' Rg (k Rg RW* - k' Rg R* g ). 
Other terms of the system account for 

• linear diffusion of W (W*) with rate D (£>*): -Dd 2 xx W (—D*d xx W*), 

• degradation of W (W*) with rate 7 (7*): -jW (-7*!^*), 

• internalization (endocytosis) of R (R*,R* g ) with rate a (a*) : —aR (-a*R*,-a*R* g ), 

• secretion of W with rate s from the source localised at x = G J: sS (5 denotes the Dirac 
Delta), 

• production of receptors: T. 

For simplicity we assume that G and F are given, strictly positive constants. 
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1.2 Nondimensionalisation 



To analyze system [HKCSj.lD we rewrite it in a nondimensional form: 

d t u\ - dl x ui = -(bi + ci + u 3 )ui + c 2 u 2 + c 4 w 4 + p%5, (t, x) G loo (2a) 

d t u 2 - dd 2 xx u 2 = -{b 2 + c 2 + c 3 u 3 )u 2 + c x u\ + c 5 w 5 , (t, x) G 1^ (2b) 

<9 t W 3 = -(63 + Mi + C 3 W 2 )W3 + C4M4 + C5M5 + P3, (*, X) E loo (2c) 

d t u A = -(64 + c 4 )m 4 + U1U3, (t, x) G loo (2d) 

= -(65 + C 5 )«5 + C 3 M 2 «3, (*, x) G /oo (2e) 



with boundary and initial conditions 

d x ui = d x u 2 = 0, 
u(0, •) = u , 

where 

u(t,x) = { Ul ,u 2 ,u 3 ,u 4 ,u 5 ){t,x) = K{W,W*,R,R*,R* g ){Tt,Lx), 
Uq(x) = (u w ,u 2 o,u 3 q,U4o,u 50 )(x) = K(W Q ,WQ,Ro,R* ,R* g0 )(Lx), 
T = L 2 /D, K = Tk R , d = D*/D, 
b = (Tj,T r y* ,Ta,Ta* , Ten*), 
c = (TkG, Tk>, k Rg /k Rl Tk' w Tk' Rg ), 
p = (KTs,0,KTT, 0,0). 



(t,x)e(dl) 

x G I 



1.3 Overview. 

The aim of this paper is to establish well-posedness of ^ and the existence of a unique steady 
state. 

During the analysis of the system ^ we encounter the following difficulties: 



absence of diffusion in equations, (2c), (2d), (2e) so that there is no smoothing effect for 
u 3 ,u 4 ,u 5 . 



singular source term in (2a), 



• nonsymmetric zero order part of the operator for the stationary problem. 

We first solve the stationary problem for ^ by using Schauder's fixed point theorem. The key 
observation is that the linear operator which appears in the definition of T n (v) (see proof of Theorem 
1), has a diagonally dominant structure. This leads us to analyze the problem in an Li(I) setting 
rather than L 2 (I). To prove uniqueness we consider the system which is satisfied by the difference of 
two possible solutions and after algebraic manipulations show that it also has a diagonally dominant 
structure. 



To remove the singularity piS from (2a) we change variables z = u — u*, where u* is the steady state 
to ([2]). Then local well-posedness in the space of continuous functions of the system for z follows 
from the classical perturbation theory for sectorial operators. To prove global existence we notice 
that the quasipositivity of the vector field appearing on the right hand side of ^ guarantees that 
the semiflow generated by ^ preserves the positive cone. Then using compensation effects it is easy 
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to show that U3, M4, w 5 G L^O, T max ; C(I)) and u±, u 2 G £oo(0, T max ; L\(I)). Finally using smoothing 
effects of the heat semigroup we prove that Ux, u 2 G Loo(0, T max ; C(J)), from which we conclude that 
system Q is globally well-posed and has bounded trajectories. 

Before stating the results precisely we introduce the notation and function spaces in which we will 
analyze the system 



2 Notation and function spaces 

For x, y G M, x, y G lR n we denote 

x V y — max{x, y}, x Ay = mm{x, y}, x + = x V 0, = (— x) V 0, sgn{x) 

n 

x = max{xj : 1 < i < n}, x = min{xj : 1 < i < n}, < x, y >= 



\x\/x , x 7^ 
,x = 



If V is a vector space we denote by V n the n-th product power of V. If (V, >) is a partially ordered 
vector space we denote its positive cone by V + := {v E V \ v > 0}. If V, W are normed spaces 
we denote by £(V, W) the space of linear, bounded operators from V to W with the usual uniform 
convergence topology. We put C(V) = C(V, V). 

We make standard convention that C denotes a positive constant which depends on various param- 
eters which are specified explicitly in the text. 

To analyze the problem we will use the following Banach spaces 

X = X = C(I), X, = C 2 N (I) = {u:ue C 2 (7), u'(-l) = = 0}, 
Xx/2 = C°'\I) = Lip(I) = W^I), 

Y = Y = Lx(I), Yx = Wl N (I) = {u:ue W?(I), u'(-l) = u'(l) = 0}, 
Y x/2 = Wl{I) = AC{I). 

Notice that due to the imbedding W?(I) C C 1 (J) the boundary conditions in the definition of Y\ are 
meaningful. 



3 Results 



From now on we assume that 



d, b > 0, c, p > 0, u G X] 



We start with the analysis of the stationary problem and prove that there exists unique nonnegative 
steady state. Observe that due to the absence of diffusion in (2c),(2d),(2e) the stationary problem 



reduces to the system (19) (see below) of two semilinear elliptic equations for u\ and u\. 
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Theorem 1. System (|2| possesses a unique nonnegative steady state 
u* G X1/2 x X\ x Xf, 2 suc/t i/iai 

U3 = F(m*,u 2 ), b 4 ul = kxulH^ul), b 5 u* 5 = k 2 u* 2 H(ul,u* 2 ), 

where 

h = 64/ (64 + c 4 ), ^2 = C365/ (65 + c 5 ), i?(:ci, ar 2 ) = Ps/ (^1 + fc 2^2 + 63 
and (u*,^) a solution of the following boundary value problem 
—it*" + (61 + ci + kxH{u\, u* 2 ))u[ - c 2 u* 2 = pi8, 



-du 2 " - ciu\ + (6 2 + c 2 + k 2 H(u\, u* 2 ))u 2 = 0, 



Ui = «, = U, 



x e I 
x e I, 
x g 9/. 



z.e. /or ever?/ <£> G Xt/ 2 



y [wtV + ((&i + ci + k%H(ul, u 2 ))u\ - 02^2)^] = si<p(0) 



and (5b) zs satisfied in the classical sense. Moreover 

u\ +pi\x\/2 G Xx. 

A typical shape of the steady state is to be found in Figure 1. 



(3a) 

(4) 

(5a) 
(5b) 
(5c) 



(6) 






Figure 1: Graph of u* (normalised to 1) computed for the following values of parameters: 
b = [100,10,10,10,10], c = [10,10,1,10,10], p = [100,0,100,0,0], d = 1/10. First row - 
M */ll u ilU,^/||^||x,^/||«3lU, second row - m^/Hm^Hx^s/IImsIIx- Notice the surprising difference 
in behavior of u\ and u* b near x = (see Remark ([I]) for explanation). 

In the following remark we analyze the behavior of the stationary solution near the source of mor- 
phogen. 

Remark 1 (u* near x = 0). Observe that as u* is unique it must be even. Indeed otherwise u*(— x) 
would be a second solution as the system ^ is invariant under the transformation x — > —x. Thus 
using ([6} 

K)'(0 + ) = - Pl /2 < 0, (7a) 
(«3)'(0) = 0, (7b) 
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where (w^)'(0 + ) denotes the right-sided derivative of u* at x = 0. 
Using (|3| and (|7| we compute directly 



I *\'/n+ 



(0 + ) = ^[^K(0),^(0))] 2 >0, (8a) 

K)'(o+) = -^[ffK(o),«;(o))] 2 K(o) + 6 3 A 2 ) < o, (8b) 

K)'(0 + ) = ^^[^(^(0),^(0))] 2 «5(0) > 0. (8c) 



In particular from (8b), ( |8c| ) and the fact that u* is even we infer that 144(0) (1*5(0)) is a strict local 
maximum (minimum), which explains the difference near x = in u\ (steep spike) and u* b (depletion 
effect) as observed in Figure ([TJ. 

Next we turn our attention to the evolution problem and establish its well-posedness and the uniform 
boundedness of trajectories in X 5 . 

Theorem 2. System (J2J) possesses a unique, global in time, nonnegative solution 

I4i G C([0, 00); X) n ^((0, 00); X) n C((0, oo); X 1/2 ) (9a) 
t4 2 e C([0, 00); X) n ^((0, 00); X) n C((0, oo); X x ) (9b) 
us, 114, U5 G ^([0, 00); X) (9c) 

stic/i t/iat /or every y> G X1/2, £ G (0, 00) 

J d t uiip + D J d x uid x (p = J [-(bi + ci + u 3 )ui + c 2 u 2 + c 4 m 4 ]<^ + Pi^(0) 

and other equations are satisfied in the sense of X . Moreover u G -^^(OjOojX 5 ) and i/ie following 
estimates hold 

5 5 

U M ^ e ~~ E ^° +^(1 - e-^/b, (10a) 

i=3 i=3 

£ lh(t)lk<e- fet ^ |Ko||y+Pi(l-e-^)/b. (10b) 

ie{i, 2,4,5} te{i, 2,4,5} 



We conclude with a remark concerning the discussion about the asymptotic behavior. 

Remark 2 (Asymptotics) . For the case of morphogen Dpp acting in the imaginal wing disc of the 
fruit fly without the presence of glypicans, it is proved in jB] that the morphogen gradient (i.e. steady 
state of the appropriate evolution system) is globally exponentially stable. It is expected that an 
analogous result should hold for [HKCSj.lD, though we are not able to prove even the local stability 
of the steady state. However it may also be the case that the presence of glypicans within a certain 
range of parameters has a destabilising effect on the equilibrium. 
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4 Lemmas 



In this section we collect lemmas which are used in the proofs of the results. The first lemma concerns 
well-posedness of the abstract ODE's in Banach spaces, while the second states that realisations of 
one dimensional Laplace operator in the chosen Banach spaces are sectorial. Since these lemmas are 
well known we state them only to make the article more self-contained. For the proofs we refer the 
interested reader to [IT] , chapter 6.1 and [12], chapter 3.1. 

Lemma 1. Assume that 

1. V\ C V are arbitrary, densely imbedded Banach spaces. 

2. A : V D V± — > V is a sectorial operator. 

3. F : V — > V is Lipschitz on bounded subsets. 
Then for given t> G V, the Cauchy problem 

v'-Av = F{v), t>0 (11a) 

v(0)=v (lib) 

has a unique, strong, maximal solution 

v G C([0, T max ); V) n ^((0, T max ); V) fl C((0, T max ); V 1 ). 

The following Duhamel formula holds 

v(t) = e At v + [ e A ^F(v(s))ds, t e [0,T max ). 
Jo 

Moreover T max is characterised by the following blow-up condition: 

T m ax < oo implies limsup||t> (t)\\v = oo. (12) 

" 7 ^ max 

In particular if there exists C > such that for t G [0, T max ) 

\\F{v{t))\\v<C{\\v{t)\\v + l), (13) 

then T max oo. 

For Z G {X, Y} we define the Z-realisation of the Laplace operator with Neumann boundary condi- 
tions on J: 

A z : Z D Z x -)> Z, A z u = u", u G Z x 

Lemma 2. Az is a sectorial, densely defined operator with compact resolvent. It generates an 
analytic, strongly continuous semigroup e tAz and for t > the following estimates hold 

\\e tA *\\£ {z) < 1, He^Hrcv^ < C7(l A t) _1/a - 
Moreover (Ax,e tAx ) is a restriction of (Ay,e tAY ) to X i.e. 

A x u = A Y u, ueXx, e tAx u = e tAy u, (t, u) G [0, oo) x X. 
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The third lemma concerns solvability of linear elliptic systems with diagonally dominant zero order 
term. It is crucial in the proofs of existence and uniqueness of the steady state of the system 



Lemma 3. Assume that for i,j = 1, 2, dj > 0, G X + and 

«n — «2i > 0, 

«22 - «12 > 0. 



(14a) 
(14b) 



Define operators 



r2 _^ y2^ 

G:Y 2 dY? -> Y'\ G = (d x A Y ) x (d 2 A Y ) + M. 



M :Y ->■ Y , Mu = (-on«i + a 12 u 2 , a 2 iu 1 - a 22 u 2 ) 

r2 -t^2 . ta2 



Tnen G is a sectorial, densely defined operator with a compact resolvent R(X, G) — (A — G) 1 and 
the following hold 



(0,oo) C p(G) and \\R(\,G)\\c(y*) < 1/A, 
ll^(A,G)|| £(y2iyi2) <C(l + l/A), 



R(X,G) preserves Fj?, 
where A > and C depends only on dj, ||ay||x • 



(15a) 
(15b) 
(15c) 



Proof. To prove that G is sectorial and has a compact resolvent notice that it is a perturbation of the 
operator (diAy) x (d 2 Ay) having these two properties by a bounded operator M G £(F 2 ). From the 
compactness of the resolvent of G we get that the spectrum a{G) only contains eigenvalues (Theorem 
6.29 0). 

In the rest of the proof we will use the following observation. Let 7 : R — > R be a function such that 



#7(x) > 0, I7I < 1 . Then using (14) we obtain the following pointwise inequality 

Mu, (7(ni),7(n 2 )) 



—0111117(^1) + ai2n 2 7(ni) + a 2 ini7(n 2 ) — a 22 u 2 ^y(u 2 ) (16a) 
= -«i7(wi)(on - 0217(^2)7(^1)) - u 2 j(u 2 )(a 22 - a i2 7(ni)7(n 2 )) < 0. (16b) 
Choose A > 0, f eY 2 ,u eYf such that 

/ = (A - G)u. 



(17) 



To proove ( 15a ) we estimate 



\U ya 



8=1 



Mm, (sgn(ui), sgn{u 2 ))) > A||n||y2 



where we used (16) with 7 = sgn and the following Kato's inequality (see Lemma 2 in [2]) 

- J v"sgn{v) >0, ueKi. (1? 



To prove (15b) observe that from (17), (15a) we have 

\\Gu\\y,< ||/|| y3 +A|H|y a <2||/|| ya , 
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whence 



u\\ Y ? < C{\\{A Y ui, A y u 2 )\\y* + ||u||ya) < C[{dx A d 2 )~ l \\Gu - Mu\\ Y 2 + ||/||WA] 

< C[(d x A d 2 )-\\\Gu\\ Y . + ||M|U {y2) H| y2 ) + \\f\\ Y 2/\] 

< C[(di A d 2 y\2 + ||M|U (y2) /A) + l/A]||/|| ya 
<C(l + l/A)||/|| y2 , 



Finally to prove (15c) assume that / G Y+. Let 4>(x) = (sgn(x) — l)/2. Then —1 < <p(x) < and 
x(j)(x) = X-. Using (16) with 7 = and (18) we obtain 

0> f (/,(0K),0M)) = A^ f Ui<j>{ui) ~l^ di [ <sgn( Ui ) - f (Mu,^^),^)) 
j i 1 i=1 J i i=1 J i J i 



> X\\u. 



\Y 2 , 



whence u > 0. 



□ 



5 Proof of Theorem 1 



We divide the proof into two parts. To prove existence of a solution of (|5]) we first approximate the 
singular source term p\5 by more regular functions h n G Y + . Using Schauder's fixed point theorem 
we prove solvability of the approximated problem. Finally using compactness methods we show that 
the approximated solutions converge to a solution of ([5]). In the proof of uniqueness we show that 
the difference of any two possible steady states belongs to the kernel of a certain operator A — G, 
where A > and G satisfies assumptions of Lemma d3| . 



5.1 Existence 

Choose a sequence h n G Y + such that h n ^* S in A^([— 1, 1]) - the space of signed Radon measures. 
For v G X+ consider the following problem 

—u" + (hi + ci + kiH{v 1: v 2 ))ui - c 2 u 2 = pih n , x E I (19a) 

—du'2 - CiU\ + (b 2 + c 2 + k 2 H(v 1 ,v 2 ))u 2 = 0, x G J, (19b) 

u[ = u ' 2 = 0, x G dl, (19c) 



where H is defined in Q. Using notation introduced in Lemma ^ system (19) is equivalent to 

(A - G)(u!, u 2 ) = (pih n , 0), 



where 



A = b, 

di — 1, d 2 — d, 

an = h -b + ci + kiH(v l ,v 2 ), a 12 = c 2 , 

021 = ci, a 22 = b 2 - b + c 2 + k 2 H(vi,v 2 ) 
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Oserve that condition (14) holds, thus using Lemma ^ we obtain that (19) has a unique solution 
(ui,u 2 ) G Y 2 + and there exists G\ which do not depend on (vi,V2), (u\,u 2 ), h n such that 



Using the compact imbedding 



||(Wl,W 2 )||y2 < Ci\\h n \\ Y . 



Y 1 CC X 



and (20) we obtain that there exists C 2 such that 

IK^i, m 2 )||x 2 < C 2 ||(ni, 1*2)11^ < CiC 2 \\h n \\ Y . 



(20) 



(21) 



(22) 



Define 



Vn = {(V1,V 2 ) G X 2 + : Wiv^V^Wx* < d^H^lly}, 

T n '-V n -> V n , T n (vi,v 2 ) = (ui,u 2 ), 



where (ui,u 2 ) is the solution of (19). Observe that V n is a closed and convex subset of a Banach 



space X 2 and T n is well defined, compact (because of (21), (22)) and continuous (because if is a 
globally Lipschitz continuous function on IR+). Hence by Schauder's theorem T n has a fixed point 

K,iX, 2 ) e v n . 



Since (h n )™ =1 is bounded in Y we get, by (22), that (^1,^2)^=1 * s bounded in Y 2 . From the 
imbeddings Y\ C Xi/ 2 CC X there exist {u\,u"^) G X\/ 2 and a subsequence (u^ ^m^ j)^ such that 
for % = 1,2 



[u 



<)', in 1^(1) 



Fix <p G X 1/2 , then since T nk (u* nk>1 , < fcj2 ) 



K fc ,i'< fc , 2 ) we have: 



y« fc ,i) v + + °i + k i H « k ,n < fe ,2)x fc ,i - c 2< fe , 2 ]^ = s i y v-^ ; 

rf y« fc , 2 ) V + [- c i< fc ,i + O2 + c 2 + k 2 H(u* nk>1 , u* nki2 ))u* nkj2 ](p = 0, (23b) 



Using again the fact that if is globally Lipschitz continuous on Mi. we can pass in (23) with n& — > 00 



and obtain that (m^,m 2 ) is a solution of (19) 



5.2 Uniqueness 

Assume that (ui,u 2 ), (vi,v 2 ) are two solutions of Noting zi = U\ — Vi for % = 1,2 we have: 

— z'l + (61 + ci)zi - c 2 z 2 + ki(H(ui,u 2 )ui - H{viiV 2 )vx) = 
-dz 2 - c 1 z l + (b 2 + c 2 )z 2 + k 2 (H(ui, u 2 )u 2 - H(v l7 v 2 )v 2 ) = 

Define 

D = (kiui + k 2 u 2 + b 3 )(k x v 1 + k 2 v 2 + 63) 
Wi = (Ui + Vi)/2, i = 1,2 
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and compute 

mv 2 - u 2 Vi = zi(u 2 + v 2 )/2 - z 2 {u\ + fi)/2 = zi%v 2 - z 2 w\ 

H(u h u 2 )ui - H(vx,v 2 )vi = p 3 (- — — - — ) = ^{k 2 (u x v 2 - u 2 v{) + 6 3 2i) 

\k\Ui + k 2 u 2 + b 3 k 1 v 1 + k 2 v 2 + b 3 / D 



j^((k 2 w 2 + 63)21 - k 2 w x z 2 ) 



H{ux,u 2 )u 2 - H(vx,v 2 )v 2 = p 3 (- — - — — j = ^{-ki{uiv 2 - u 2 v x ) + b 3 z 2 ) 

\faut + k 2 u 2 + 3 kiVi + k 2 v 2 + 63/ D 



^{-k x w 2 zx + {k\W\ + 63)2:2). 



Thus 



// , /, . hp 3 b 3 . . hk 2 p 3 w 2 hk 2 p 3 wi 
~H + (h + — ^— + ci + ^ > Zl ~ (° 2 + d ' Z2 = 

, „ , . hk 2 p 3 w 2 k 2 p 3 b 3 kik 2 p 3 wi 
-dz 2 - (ci H — )z 1 + (6 2 H — h c 2 H — )2 2 = 

Hence, using the notation introduced in Lemma ([3]), (21,2:2) belongs to the kernel of the operator 
h — G where 

di = 1, d 2 — d 

, , , &iP3&3 , , hk 2 p 3 w 2 k 1 k 2 p 3 w 1 
a 1 i = 6i-bH — hciH — , a 12 = c 2 H — , 

hk 2 p 3 w 2 k 2 p 3 b 3 k x k 2 p 3 wx 
a>2\ = c\ H ^ , a 2 2 = 6 2 -bH — h c 2 H — . 



Since nonnegativity of Wi, w 2 ensures that assumption (14) is fulfilled we infer that 21 = 2:2 = which 
finishes the proof. 

5.3 Regularity of u\ outside x = 0. 



Observe that E = —p\\x\/2 satisfies —E" = p\5 in the sense of distributions. Owing to (5a) v = u\— E 
solves the following boundary value problem 

-v" = f, xel 

v = —E', xedl 
with / = C2W2 — (^1 + Ci + k\H{u\, u 2 ))u*. Since / 6 X then g follows. 



6 Proof of Theorem 2 

Using the theory of analytic semigroups we first establish the local well-posedness of (J2J). Using 
quasipositivity of the right hand side of (|2j) we next prove that the generated semiflow preserves 
nonnegativity of initial conditions. Then using a compensation effect we derive Loo(0, oo, X) estimate 
for m 3 ,m 4 ,m 5 and L^^O, oo,Y) estimate for Ux,u 2 . Finally thanks to the regularising properties of 
the semigroup e Ayt we bootstrap the estimate to u 6 Loo(0, oo,X 5 ). 
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6.1 Local existence 



We rewrite system ^ in the new variables z = u — u*, where u* is the unique steady state of ([2]), 
and put it into the semigroup framework: 



t! - A x z = f(z), 

z(0) = z = u - u* 



t > 



(24a) 
(24b) 



where 



A x 
f 

/ 2 (z 
/ 3 (z 

; 4 (z 

/ 5 (z 



Ax x (dA x ) x 3 

(/,../2../,./;../^:.V 5 ^.V^ 

-(&i + c\)zi - (ziz 3 + u\z 3 + U* 3 Zi) + c 2 z 2 + c A z 4 

-{h + c 2 )z 2 - c 3 (z 2 z 3 + u* 3 z 2 + u* 2 z 3 ) + c x z x + c 5 z 5 

-b 3 Z 3 - ( y Z 1 Z 3 + U[Z 3 + U*^) - C 3 (Z 2 Z 3 + U* 3 Z 2 + U* 2 Z 3 ) + C424 + C5Z5 
-(64 + C 4 )^4 + {Z\Z 3 + U[Z 3 + W324) 
-(&5 + C5)25 + C 3 (z 2 Z 3 + U* 3 Z 2 + U* 2 Z 3 ). 



Observe that Ax generates an analytic, strongly continuous semigroup in X 5 
e* Ax = e tAx x e tdAx x (Id) 3 . Moreover f is Lipschitz continuous on bounded subsets of X 5 . Using 



Lemma Q we obtain that (24) possesses a unique solution defined on a maximal time interval 
[0, T max ) with the following regularity: 



z 1 ,z 2 eC([0,T max );X)nC 1 ((0,T n 
Z3, Z4, Z5 



X)nC((0,T max );X 1 ) 



Setting u = z + u* it is obvious that u is the unique solution to (J2|. 



6.2 Positive cone invar iance 

Consider the following system 



dtvi - d 2 xx vx = 


-(61 + ci + (f 3 )+)ui + c 2 (v 2 ) + + c A (v A )+ + p ± 5, 


(t,x) 


e ^oo 


(25a) 


d t v 2 - dd 2 xx v 2 = 


-{b 2 + c 2 + c 3 (v 3 ) + )v 2 + ci(fi)+ + c 5 (v 5 ) + , 


(t,x) 


e loo 


(25b) 


d t v 3 = 


-(63 + («l)+ + C 3 (f2)+)^3 + c 4 (f 4 )+ + c 5 (w 5 ) + + p 3 , 


(t,x) 


e loo 


(25c) 


d t v A = 


-(64 + C 4 )^4 + (Vx) + (vz) + , 


(t,x) 




(25d) 


d t v 5 = 


-(h + c 5 )v 5 + c 3 (v 2 )+(v 3 ) + , 


(t,x) 


e /oo 


(25e) 



with boundary and initial conditions 



d x vi = d x v 2 = 0, (t, x) E [8I) C 

v(0, •) = u , x e I 



Reasoning as in the previous section, the system (25) possesses a unique maximally defined solution 



v(i) on [0,T^ ax ) in X 5 . We will now prove that v(t) > for t G [0,1^). Multiplying (25) by 
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— (v)_ and adding equations we obtain: 



d/dtJ2\\(vi)-\\l + WdM-Wl + 2d\\d x (v 2 u\l < o 

D(«<(o))-||2 = o. 



1=1 



Thus v(£) > for t G [0,T^ ax .). Then v + = v and it readily follows from (25) that v solves (|2]) on 
[°> r moJ- Consequently, v = u on [0,T^ ax ) and T^ ax . < T max . Finally observe that if T' max < oo 



then, by (12), lim sup^y/ ||u(t)||x5 = limsup t _j. T / ||v(t)||x s = oo thus T' max = T max and u(t) > 
on [0,T max ). 



6.3 w 3 , U4, w 5 G 

-^00(0? T max , 



Adding equations (2c),(2d),(2e) and using nonnegativity of u, we obtain for every x £ I 

5 5 



dt^Ui + b^m <p 3 , te[0,T n 



i=3 i=3 



Thus 



0<Y, u i< e ~~ Yl Ul ° + M 1 - e ~~') A t e [0, T n 



(26) 



i=3 



i=3 



6.4 T max 00 



Observe that, due to (26), f(z(t)) satisfies (13) with V = X 5 hence T max = 00 by Lemma Q 



6.5 Ui, U 2 G ^oo(0, 00; Y) 



After integrating equations (2a), (2b), (2d), (2e) over the set / and adding them together we ob- 
tain 



d_ 



\\ui\\ Y )+b ^2 \\ u i\w<pi, 



ie{l, 2,4,5} 



ie{l,2,4,5} 



Thus 



X] IK(t)||y<e^ ]T ||«io||y+Pi(l-e-^)/6. 

ie{l, 2,4,5} ie{l, 2,4,5} 



(27) 
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6.6 U\, U2 6 ^oo(0, oo; X) 



From (26), rt27| we obtain that /i(z) +Zj £ £oo(0, oo; Y). Using the Duhamel formula and estimates 



from Lemma (]2|) we get 

\\zi(t)\\ x < e-'We^WaxMoWx + [ e~ s \\e sAY \\ c( y,x) WfMt - s)) + z x (t - s)\\ Y ds 

Jo 

POO 

< \\z 10 \\x + C (1 A sy^e-dsWhiz) + zi\\ Loo <y) < IkiolU + C||A(z) + Zl \\ Loo(Yh 
Jo 

whence U\ G Loo(0, oo; X). A similar argument gives Ui G Loo(0, oo; X) and completes the proof. 
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